Homework 12
Please Add Your Name Here
2024-03-05
Instructions
Answer the following questions and complete the exercises in
RMarkdown. Please embed all of your code and push your final work to
your repository. Your final lab report should be organized, clean, and
run free from errors. Remember, you must remove the # for
the included code chunks to run. Be sure to add your name to the author
header above. For any included plots, make sure they are clearly
labeled. You are free to use any plot type that you feel best
communicates the results of your analysis.
Make sure to use the formatting conventions of RMarkdown to make your report neat and clean!
Load the libraries
library(tidyverse)
library(janitor)
library(ggmap)Load the Data
We will use two separate data sets for this homework.
register_stadiamaps("4b286b86-a5e8-4791-a641-ac6a0aede71a", write = FALSE)- The first data set represent sightings of grizzly bears (Ursos arctos) in Alaska.
bears <- readr::read_csv("data/bear-sightings.csv") %>% clean_names## Rows: 494 Columns: 3
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## dbl (3): bear.id, longitude, latitude
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.- The second data set is from Brandell, Ellen E (2021), Serological dataset and R code for: Patterns and processes of pathogen exposure in gray wolves across North America, Dryad, Dataset.
wolves <- readr::read_csv("data/wolves_data/wolves_dataset.csv") %>% clean_names## Rows: 1986 Columns: 23
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (4): pop, age.cat, sex, color
## dbl (19): year, lat, long, habitat, human, pop.density, pack.size, standard....
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.- Load the
grizzlydata and evaluate its structure.
summary(bears)## bear_id longitude latitude
## Min. : 7 Min. :-166.2 Min. :55.02
## 1st Qu.:2569 1st Qu.:-154.2 1st Qu.:58.13
## Median :4822 Median :-151.0 Median :60.97
## Mean :4935 Mean :-149.1 Mean :61.41
## 3rd Qu.:7387 3rd Qu.:-145.6 3rd Qu.:64.13
## Max. :9996 Max. :-131.3 Max. :70.37- Use the range of the latitude and longitude to build an appropriate bounding box for your map.
long <- c(-166.2,-131.3)
lat <- c(55.02,70.37)
bbox <- make_bbox(long, lat, f = 0.03)- Load a map from
stamenin a terrain style projection and display the map.
bearsmap <- get_stadiamap(bbox, maptype = "stamen_terrain", zoom=7)## ℹ © Stadia Maps © Stamen Design © OpenMapTiles © OpenStreetMap contributors.## ℹ 196 tiles needed, this may take a while (try a smaller zoom?)- Build a final map that overlays the recorded observations of grizzly bears in Alaska.
ggmap(bearsmap) +
geom_point(data=bears, aes(longitude, latitude), size=0.4)
Let’s switch to the wolves data. Brandell, Ellen E (2021), Serological dataset and R code for: Patterns and processes of pathogen exposure in gray wolves across North America, Dryad, Dataset.
- Load the data and evaluate its structure.
summary(wolves)## pop year age_cat sex
## Length:1986 Min. :1992 Length:1986 Length:1986
## Class :character 1st Qu.:2006 Class :character Class :character
## Mode :character Median :2011 Mode :character Mode :character
## Mean :2010
## 3rd Qu.:2016
## Max. :2019
##
## color lat long habitat
## Length:1986 Min. :33.89 Min. :-157.84 Min. : 254.1
## Class :character 1st Qu.:44.60 1st Qu.:-123.73 1st Qu.:10375.2
## Mode :character Median :46.83 Median :-110.99 Median :11211.3
## Mean :50.43 Mean :-116.86 Mean :12797.4
## 3rd Qu.:57.89 3rd Qu.:-110.55 3rd Qu.:11860.8
## Max. :80.50 Max. : -82.42 Max. :34676.6
##
## human pop_density pack_size standard_habitat
## Min. : 0.02 Min. : 3.74 Min. :3.55 Min. :-1.63390
## 1st Qu.: 80.60 1st Qu.: 7.40 1st Qu.:5.62 1st Qu.:-0.30620
## Median :2787.67 Median :11.63 Median :6.37 Median :-0.19650
## Mean :2335.38 Mean :14.91 Mean :6.47 Mean : 0.01158
## 3rd Qu.:3973.47 3rd Qu.:25.32 3rd Qu.:8.25 3rd Qu.:-0.11130
## Max. :6228.64 Max. :33.96 Max. :9.56 Max. : 2.88180
##
## standard_human standard_pop standard_packsize standard_latitude
## Min. :-0.9834 Min. :-1.13460 Min. :-1.7585 Min. :-1.805900
## 1st Qu.:-0.9444 1st Qu.:-0.74630 1st Qu.:-0.5418 1st Qu.:-0.636900
## Median : 0.3648 Median :-0.29760 Median :-0.1009 Median :-0.392600
## Mean : 0.1461 Mean : 0.05084 Mean :-0.0422 Mean :-0.000006
## 3rd Qu.: 0.9383 3rd Qu.: 1.15480 3rd Qu.: 1.0041 3rd Qu.: 0.814300
## Max. : 2.0290 Max. : 2.07150 Max. : 1.7742 Max. : 3.281900
##
## standard_longitude cav_binary cdv_binary cpv_binary
## Min. :-2.144100 Min. :0.0000 Min. :0.0000 Min. :0.0000
## 1st Qu.:-0.359500 1st Qu.:1.0000 1st Qu.:0.0000 1st Qu.:1.0000
## Median : 0.306900 Median :1.0000 Median :0.0000 Median :1.0000
## Mean :-0.000005 Mean :0.8529 Mean :0.2219 Mean :0.7943
## 3rd Qu.: 0.330200 3rd Qu.:1.0000 3rd Qu.:0.0000 3rd Qu.:1.0000
## Max. : 1.801500 Max. :1.0000 Max. :1.0000 Max. :1.0000
## NA's :321 NA's :21 NA's :7
## chv_binary neo_binary toxo_binary
## Min. :0.0000 Min. :0.0000 Min. :0.0000
## 1st Qu.:1.0000 1st Qu.:0.0000 1st Qu.:0.0000
## Median :1.0000 Median :0.0000 Median :0.0000
## Mean :0.8018 Mean :0.2804 Mean :0.4832
## 3rd Qu.:1.0000 3rd Qu.:1.0000 3rd Qu.:1.0000
## Max. :1.0000 Max. :1.0000 Max. :1.0000
## NA's :548 NA's :538 NA's :827- How many distinct wolf populations are included in this study? Mae a new object that restricts the data to the wolf populations in the lower 48 US states.
wolves %>%
summarise(distinct= n_distinct(pop))## # A tibble: 1 × 1
## distinct
## <int>
## 1 17- Use the range of the latitude and longitude to build an appropriate bounding box for your map.
wolves %>%
select(lat, long) %>%
summary()## lat long
## Min. :33.89 Min. :-157.84
## 1st Qu.:44.60 1st Qu.:-123.73
## Median :46.83 Median :-110.99
## Mean :50.43 Mean :-116.86
## 3rd Qu.:57.89 3rd Qu.:-110.55
## Max. :80.50 Max. : -82.42lat <- c(33.89, 80.50)
long <- c(-157.84, -82.42)
bbox2 <- make_bbox(long, lat, f = 0.03)- Load a map from
stamenin aterrain-linesprojection and display the map.
wolvesmap <- get_stadiamap(bbox2, maptype = "stamen_terrain", zoom=7)## ℹ © Stadia Maps © Stamen Design © OpenMapTiles © OpenStreetMap contributors.## ℹ 1218 tiles needed, this may take a while (try a smaller zoom?)- Build a final map that overlays the recorded observations of wolves in the lower 48 states.
ggmap(wolvesmap)+
geom_point(data=wolves,aes(x=long, y=lat),size=0.4)
- Use the map from #9 above, but add some aesthetics. Try to
fillandcolorby population.
ggmap(wolvesmap)+
geom_point(data=wolves,aes(x=long, y=lat, color=pop),size=0.4)
Push your final code to GitHub!
Please be sure that you check the keep md file in the
knit preferences.